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We study the tunneling density of states (TDOS) for a junction of three Tomonaga-Luttinger liquid 
wires. We show that there are fixed points which allow for the enhancement of the TDOS, which is 
unusual for Luttinger liquids. The distance from the junction over which this enhancement occurs is 
of the order of x — v/(2u), where v is the plasmon velocity and lo is the bias frequency. Beyond this 
distance, the TDOS crosses over to the standard bulk value independent of the fixed point describing 
the junction. This finite range of distances opens up the possibility of experimentally probing the 
enhancement in each wire individually. 
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PACS numbers: 71.10.Pm, 71.27.+a, 73.40.Gk 



Introduction.- Junctions of multiple quantum wires 
(QWs) have attracted considerable attention in recent 
years since they form the basic building blocks of quan- 
tum circuitry. Experimental realizations of three-wire 
Y-junctions of carbon nanotubes have also given the 
field a strong impetus pj. Earlier studies of junctions 
of Tomonaga-Luttinger liquids (LL) were mostly focused 
on searching for various low-energy fixed points (FPs) 
@- H, IH H, W, S B @] an d the corresponding conductance 
matrices; this includes the recent development involving 
the inclusion of a finite superconducting gap at the junc- 



tion 



Here we will focus on the local single-particle tunneling 
density of states (TDOS); this can be found by measur- 
ing the differential tunneling conductance of a scanning 
tunneling microscope (STM) tip at a finite bias. The 
differential conductance provides a direct measure of the 
TDOS of the electrons at an energy given by the bias 
voltage [13], provided the density of states for the STM is 
energy independent. For a two-wire LL junction, the cur- 
rent measured by the STM tip varies as a power of the 
bias [13] ; the power depends on the Luttinger parameter 
g and the FP to which the junction has been tuned. 

Earlier studies of the TDOS in a LL system with 14] 
and without [l5[ impurities revealed that the TDOS van- 
ishes as a power law in the zero bias limit [l9j . This is 
popularly considered to be a hallmark of a LL system. 
The TDOS for a LL wire with an impurity was studied in 
Ref. [HI where the tunneling STM current close to the 
impurity was shown to be power law suppressed in the 
zero bias limit. An enhancement of the TDOS was found 
at a single LL -superconductor [lfij] junction and was ex- 
plained in terms of the proximity effect. An enhancement 
of the spectral weight was also found at a junction of 
multiple LLs tuned to a fermionic FP (linear boundary 
condition (BC) between fermion fields at the junction) 
Q. All the earlier studies of the TDOS were focused on 
either a multiple wire junction tuned to a fermionic FP or 
a single LL -superconductor junction, but not a junction 



of multiple LLs tuned to bosonic FPs. Bosonic FPs refer 
to linear BCs connecting incoming and outgoing currents 
(which are bilinears in the fermion fields) at the junction. 

Here we study the TDOS of a three-wire junction of a 
single channel QW modeled as a LL. We find that close 
to the junction, the TDOS depends on the details of the 
current splitting matrix at the junction which relates the 
incoming and outgoing currents; this matrix describes 
bosonic FPs of the system. We find that for a certain 
range of repulsive intcr-clectron interaction (g < 1) and 
certain current splitting matrices, the TDOS close to the 
junction shows an enhancement in the zero bias limit. We 
show that this is related to reflection of holes off the junc- 
tion which mimics the Andreev reflection process, even 
though there is no superconductor in the present sce- 
nario. This is in sharp contrast to the case of a two- wire 
junction where a repulsive electron-electron interaction 
always results in a suppression of the TDOS near the 
junction. This is the central result of this Letter. Far- 
away from the junction, the TDOS reduces to that of the 



bulk LL wire, p(uj) 



,( 9 + 3 - 1 -2)/2 



independent of the 



details of the junction, and shows a suppression for both 
repulsive (g < 1) and attractive (g > 1) interactions. 

Bosonization.- The electron field (taken to be spinless 
for simplicity) can be bosonized as ip(x) — (l/\/2wa) 
[F e* 1 *' 1 ' + F I e- ikFX+i ^^), where (f>o(x), F 
and (j>i(x), Fj are the outgoing and incoming chiral 
bosonic fields and corresponding Klein factors, Hf is the 
Fermi momentum, and a is a short distance cut-off. We 
model the wires as spinless LLs on a half-line (x > 0), i.e., 
all the wires are parametrized by a coordinate x running 
from to oo. The corresponding Hamiltonian is given by 



H 



dx 
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where prime (dot) stands for spatial (time) derivative, 
<f>i(x,t) = [<j>oi + <f>n)/2, Oi(x,t) = (4>oi - 0«)/2, 4>i = 
(v/g) 6' i: and 0, = (vg) The total density and current 
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= (1 ' 


h 2cos<9)/3, 


6(c) 



can be expressed in terms of the incoming and outgo- 
ing fields as p = po + pi with p Q /i = ±(l/2n) 4>' ^ I: 
and j = j - jj with j a/I = ±(v F /2ir) (f>' 0/1 . Finally, 
we need to impose a boundary condition on the fields 
at s = 0. Following Ref. 7], the incoming and out- 
going currents, and consequently the bosonic fields, are 
related at the junction by a current splitting matrix M, 
i.e., joi = J2j M ij Jij, which leads to <j) 0l = J2j M y 4>ij- 
To ensure that the matrix M represents a FP of the the- 
ory, the incoming and outgoing fields must satisfy the 
bosonic commutation relations; this restricts the matrix 
M to be orthogonal. Scale invariance imposes the same 
constraint of orthogonality on M as shown in Ref. ll| ; 
orthogonality also implies that there is no dissipation in 
the system. Current conservation at the junction implies 
that each row and column of M adds up to 1 . 

In general, for a three- wire charge-conserving junction, 
M can be parametrized by a single continuous parameter 
8, and it falls into one of two classes with (a) dct Mi = I, 
and (b) detM2 = — 1. These classes can be expressed as 



(2) 



cos# + (— )v / 3sin#)/3. This provides us with an ex- 
plicit parametrization of the two one-parameter families 
of FPs; any FP in the theory can be identified in terms 
of 8, with the FPs at 8 = and 8 = 2ir being identical. 
The det Mi = I class represents a Z$ symmetric (in the 
wire index) class of FPs, while detM2 = — I represents 
an asymmetric class of FPs. In the Mi class, 8 — ir cor- 
responds to the Dp FP, 8 = corresponds to the discon- 
nected N FP, and 8 = 2ir/3 and 4ir/3 correspond to the 
chiral cases x± > following the notation of Ref. [B| . Since 
4>o and <pi are interacting fields, we must perform a Bo- 
goliubov transformation, 4>o/i = (1/2V5){(1 + 9)4>o/i + 
(1 — g)4>i/o}> to obtain the corresponding free outgoing 
(incoming) (4>o/i) chiral fields satisfying the commuta- 
tion relations, [4>o/i{x, t), 4>o/i{ x 'i t)] = ±iirSgn(x — x'). 
Unlike the usual Bogoliubov transformation, here we 
also need to consider the effect of the junction ma- 
trix M relating the interacting incoming and outgo- 
ing fields. Following Ref. 7], we obtain a Bogoli- 
ubov transformed matrix M which relates <j>oi to <pn. 
We find that 4>oi(x) = J2j M v 4>ij(- x ) where M = 
[(l+g)I- (I-g)Mr 1 [(l + g)M- (I — <?)!]. Note that 
the M 2 class of matrices satisfy (M 2 ) 2 = I; hence M 2 = 
M 2 , and the interacting and free fields satisfy identical 
BCs at the junction. This is not true for the Mi class, 
but Mi still has the same form as Mi with the corre- 
sponding parameters a = (3g 2 — 1 + (3g 2 + 1) cos#)/ry 
and 6(c) = 2(1 — cosd + (— ) ^/5g sin 8) /rj, where r/ = 
3(1 +.9 2 + (.o 2 -I)cos0). 

The M matrix is related to the DC conductance matrix 
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FIG. 1: Schematic picture of STM tip for measuring the 
TDOS near the junction, the region of inter-wire tunnelings 
and a magnetic field B at the junction. 

iven by G = (2e 2 /h)(l - M) for Fermi liquid leads 
llj . Qualitatively, M is related to tunnelings between 
the different wires and backscatterings in each wire. The 
experimental set-up can be a junction of several edges of 
a quantum Hall system as in Ref. [7j, and M (or 8) can 
be tuned by applying gate voltages and a magnetic field 
at the junction. We note that M 2 is time-reversal invari- 
ant, but Mi is generally not and tuning it will require a 
magnetic field piercing the junction (see Fig. [TJ. 

Tunneling density of states.- We now compute the 
TDOS for adding an electron with energy hut on the i th 
wire Hal. 



2tt^|(0| 4 (:r) \ n)\ 2 5(E n - E - hw) 
= 2R< I dt (0| ipi(x,t)ipl(x,0) |0) e Mt .(3) 



Here |n) (E n ) denotes the n th eigenstate (eigenvalue) of 
the Hamiltonian in Eq. ([TJ. The Green's function in the 
i th wire is G = {^i{x,t)4{x,$)) = (ipn(x, t)^ H {x, 0)) + 
(iPo t (x, t)^ Ql {x, 0)) + e~ 2 ^ x (iPn (x, t)^ Ql {x, 0)) + e 2lk ^ x 
{fpOi(x, t)f/'j'i(^, 0)). The two non-oscillatory terms are 

(^OM)^(x,0)) = (il> Oi (x,t)i> Oi (x,0)) 







(1 + 9 2 ) 
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(4) 



The oscillatory part vanishes as L — > oo and can be 
dropped in further discussions. For the Mi class, di = a; 
for the M 2 class, di = a,b,c depending on the wire index. 

Treating the tunneling strength 7 between the i th wire 
and the STM tip perturbatively and using Eqs. (JSJUJ), the 
differential tunneling conductance evaluated to leading 
order in 7 is found to be directly proportional [l2| to the 
TDOS on the i th wire. The TDOS has the same form in 
the i->0 and x — » 00 limits and is given by 



ahT{A) 



A A- 



Aa/v e( w ), (5) 



where T(A) is the Gamma function, 9(w) is the Heav- 
iside step function, r c = a/v is the short time cut-off, 
and lo = eV/h, where e is the electronic charge and V 



3 



is the bias voltage between the STM tip and wire sys- 
tem held at a uniform potential. The cut-off frequency 
scale for the validity of the perturbation theory is given 
by ojq = [|7| _2 /( A<_1 )«]/a. For x — > 0, A is a function 
of the di which is the corresponding diagonal element of 
the appropriate M matrix. For the Mi class, A = An (5), 
while for Ma, A = A* is a function of a, b, c depending on 
the wire index i. For x —> oo, A = (g + g )/2 indepen- 
dent of the M matrix at the junction. Thus we recover 
the expression for the bulk TDOS in a LL as x — > oo 
irrespective of the details of the junction. For the Mi 
class, the power law exponent for x — > is the same on 
all the three wires due to Z% symmetry and is given by 



A = 



1 5# 2 + 1 + ( 5 2 - 1) cosfl 
3g~ g 2 + l + (g 2 -l) cos(9 ' 



(6) 



Eq. ijSJ) indicates that for g < 1 there are values of 9 
for which A < 1. This implies that there are FPs in 
the theory which show an enhancement (see Eq. (|5[l) 
of the TDOS in the zero bias limit cu — > 0. This is in 
sharp contrast to previous studies for various cases of 
normal (not superconducting) junctions of two-wire sys- 
tems which always showed a suppression of the TDOS in 
the zero bias limit for g < 1. This is our main result. 
Note that whenever g = 1, A = 1; this implies that the 
TDOS is independent of the bias for g = 1. This might 
look natural since g = 1 corresponds to free fermions 
in the wire for which the TDOS is expected to be bias 
(energy) independent. However, this is misleading; the 
BC conditions expressed in terms of the matrix M at 
the junction correspond to non-linear relations between 
the fermions on each wire in the vicinity of the junction, 
and hence represent non-trivial interaction between the 
fermions at the junction. Hence, the TDOS being energy 
independent for g = 1 for any FP represented by Mi is 
a non-trivial result by itself. To get a clear idea about 
the FPs which show an enhancement of the TDOS , we 
present contour plots of A in the g — 9 plane in Fig. [2] 
In the left plot (corresponding to the Mi class), we can 
see a dome shaped region for g < 1 which corresponds 
to FPs showing an enhancement, i.e., An < 1. It is in- 
teresting to note that this region is bounded by the two 
chiral FPs, x± at 9 = 2tt/3, 4tt/3. The D P FP at 9 = ir 
also falls in this region and shows an enhancement of the 
TDOS for 1/2 < g < 1. For the M 2 class, the power law 
exponents for the three wires are given by 



Ai 



4 + 2g 2 + {cos0- \/3sin60 (g 2 - 1) 
6g 



(J) 



A 2 and A 3 which are obtained by shifting 9 — > 9 =p 2ir/3 
in A x . For this class also, there are FPs which show an 
enhancement of the TDOS for g < 1; they correspond 
to the dome shaped region in the right plot in Fig. [5] 
In contrast to the Mi class, in this case there can be an 
enhancement in one wire and suppression in the other 
wires due to the broken Z3 symmetry. 
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FIG. 2: Ao for all wires in the Mi class (left), and Ai for wire 
1 in the M2 class (right) in the g — 9 plane. 

The observability of the different FPs and the en- 
hancement of the TDOS near the junction crucially de- 
pends on the renormalization group (RG) stability of the 
FP against various perturbations in the form of inter- 
wire electron tunnelings at the junction (see Fig. [T]). To 
understand the stability of the FP, we present the scal- 
ing dimension Sq for all possible single-electron tunneling 
events (quadratic in fermion operators) in TableslIlandHTl 
in these tables, ipi/o are related to 4>i/o by the bosoniza- 
tion formula given earlier. In Fig. [31 we plot the values 
of 1 — Sq as functions of 9 for g = 0.9. A positive value 
of 1 — Sq implies that the operator is RG relevant. 



Operator 



Scaling dimension So 



■4>20^I^lo^2ui > loi ) 3I 
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3(s a +cos 6(9^-1)4-1) 
2g(cos e+V3sinS + 2) 


3(3^+0030(9^ 
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n« + 2) 




-D+i) 
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t/ 2 +cos 6{g' 2 - 


-i)+i 




2g 




q 2 +cos 0(q' 2 - 


-D+i 



TABLE I: Table of tunneling operators for the Mi class. 



Operator 



Scaling dimension do 



V4oV>ii 
^20^2/ 



±y(2 - 2 cos 6) 
\g{2 + cos e + v^ sine) 
|g(2 + cosJj - Vising) 



i>\oi>2i,^lo^u ^T2T (2-2cos0) 
i>loi>3i^to^2i ^(2 + cost)- V3sin( 
V'aoV'i/, ^lo^si ^(2 + cosg + V3sin( 



^11^21, 1>toi>20 



2(9^ + l) + (9^-l)2cos9 
49 

,/,t ,/. ,/.t ,/, 2(9 2 + l)-(9 2 -l)(cos9-y3 S m9) 
,/,t ,/, ,/.t ,/. 2(9 2 +l)-(g 2 -l)(co5 9+y3 S in9) 

TABLE II: Table of tunneling operators for the M2 class. 



It is clear from Fig. [3] that all the FPs showing an en- 
hancement of the TDOS are unstable; however, for the 
Mi class close to 9 = 2tt/3, 4tt/3 (the x± FPs), only 
one operator is highly relevant and the rest are almost 
marginal. Hence this part of the parameter space al- 
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FIG. 3: 1 — So as functions of for the classes Mi (left) and 
M2 (right) for various tunneling operators, for g = 0.9. 

lows for a large temperature window for observing an en- 
hancement if one can experimentally suppress the most 
relevant tunneling by tuning the junction appropriately. 

To gain a better understanding of the enhancement, 
we expand the A's for g ~ 1 in the small parameter 
1 — g to obtain Aj = 1 + (1 — g)dj, where j = 0, 1, 2, 3, 
and dj are the diagonal elements of the corresponding 
M matrix. This limit corresponds to weakly interacting 
electrons in the bulk of the wires away from the junction. 
Whenever dj < and g < 1, Aj is less then unity which 
corresponds to an enhancement of the TDOS in the zero 
bias limit. But dj < corresponds to a hole current being 
reflected from the junction when an electron current is 
incident on the junction. Hence we conclude that all the 
FPs which involve reflection of a hole off the junction lead 
to an enhancement of the TDOS. As discussed earlier, an 
enhancement of the TDOS was previously observed in a 
junction of a LL wire with a superconductor (ljj ; this can 
be attributed to the proximity induced Andreev process 
at the junction which results in the reflection of a hole 
from the junction in response to an incident electron. It is 
interesting to note that for our case too, the enhancement 
is connected to holes being reflected off the junction, even 
though there is no superconductor in our model. Finally, 
note from Eqs. (HHS]) that the cross-over length scale 
beyond which the TDOS goes over to its bulk form is 
given by x = vj (2u>). For a typical bias voltage of 10 [iV 
and a Fermi velocity v ~ 10 5 m/s (typical of a two- 
dimensional electron gas), we get a cross-over length of 
about 3 fim; this is readily accessible within present day 
experimental realizations of a one-dimensional QW in a 
two-dimensional electron gas. 

Discussion.- It is important to note that the interest- 
ing prediction of an enhancement of the TDOS for g < 1 
involves unstable FPs; hence the enhancement can be 
expected to be observed in experiments at high tempera- 
tures only. If the junction is tuned to one of the FPs which 
show enhancement, a variation of the temperature from 
high to low will first show an enhancement and then a 
suppression of the TDOS as the system finally flows to 
the disconnected stable FP at low temperatures. This 
non-monotonicity observed via the STM current will be 



a hallmark of our prediction. Here, high and low temper- 
atures are defined with respect to a cross-over scale called 
wo after Eq. (0. In other Luttinger liquid systems, such 
as the one studied experimentally in Ref. [18[ , the cross- 
over scale Tb was found to be of the order of 0.5 — 3K 
(corresponding to u ~ /cbTb/S ~ 60 — 400GHz); this 
scale can be varied by tuning the tunneling strength 7. 
Thus the cross-over scale can be tuned experimentally as 
was done in Ref. [H| . Hence the temperature window in 
which the enhancement of the TDOS can be observed is 
experimentally tunable. 
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